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Abstract 

With characteristics of non-linear and time-varied, so it is 
difficult for a marine Diesel Engine to be controlled with 
traditional PID controller. An adaptive controller based on 
back-propagation (BP) neural network and Wiener model 
identifier was put forwarded to tune PID parameters for 
marine diesel engine speed control system. In the controller, 
a Wiener neural network structure was applied to identify 
Wiener model of diesel engine nonlinear model. The weights 
in the Wiener neural network are adjusted with 
backward-propagation methods, and those weights stand 
for the parameters of the Wiener model. In order to satisfy 
the different work conditions, the adaptive controller is 
improved via introducing relative error in target evaluation 
function of the BP neural network. Moreover, with the 
sensitivity function of diesel engine output with respect to 
its input obtained by the WNN identifier, the convergence 
speed of optimizing PID parameters is improved. A 
simulated test on a diesel engine demonstrated that the 
adaptive controller improved control performance over the 
traditional PID control. 
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Introduction 

Speed governing system of diesel engine is an 
important component of diesel engine electronic 
control system. At present, the traditional PID 
controller is one of the popular controllers in diesel 
engine speed governing system, since its design is 
simple and does not require detailed knowledge of the 
system dynamics [1]. However, diesel engine is a 
non-linear and time-varied system, and it is difficult to 
determine the appropriate PID parameters when 
various uncertainties and nonlinearities exist. 
Moreover, once control parameters established, it is 
difficult to be adjusted online. Numerous research 
papers focused on adaptive PID control [2], self-tuning 
PID control [3, 4], etc. In recent years, the combination 
of the PID control and neural network has become a 



new direction for intelligent controller, attracting many 
researchers [5-10]. The self-tuning PID can use the BP 
neural network to tune PID parameters online in 
accordance with change of load and conditions 
automatically. However, because partial derivative of 
diesel engine output with respect to its input cannot be 
calculated directly, convergence speed of the BP neural 
network tuning PID parameters is not fast. In order to 
increase convergence speed, it is necessary to identify 
the model of diesel engine, and obtain sensitivity 
function of diesel engine output with respect to its 
input. 

As to nonlinear system identification, Wiener models 
are widely used for modeling of processes [11-14], 
which describe a diesel engine model as a linear 
dynamic block in cascade with a nonlinear static gain. 
In order to achieve a certain performance or adaptive 
capability, neural networks were integrated into 
Wiener model to form Wiener neural network [15, 16]. 
However, a diesel engine is a time delay system, so 
time delay must be considered in Wiener neural 
network. 

In this paper, an adaptive controller is proposed, which 
includes a Wiener model identifier and a PID 
controller tuned by back-propagation (BP) neural 
network (BP-PID). The Wiener model identifier is a 
Wiener-type neural network, which is used to identify 
diesel engine model and to provide approximate 
sensitivity function of the diesel engine model online 
to update the PB-PID controller. Moreover, considering 
different work conditions, relative error is introduced 
in target evaluation function of the BP neural network. 
The adaptive controller was tested on identifying and 
controlling a diesel engine, and the experiment 
confirmed the effectiveness of the proposed method. 

The Control Model of Marine Diesel Engine 

From a control point of view, two important paths 
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have to be considered in marine diesel engine: fuel and 
air. Fig.l shows a schematic overview of the basic 
structure of a typical diesel engine control system. The 
control inputs to the fuel path are start of injection 
duration, and injection pressure. The control outputs 
to the fuel path are torque, speed and exhaust-gas 
emissions. The turbocharger dominates the air path, 
and it is controlled using a closed-loop approach 
where the measured output is the boost pressure. 



The Cause and effect diagram of a typical diesel engine 
is shown in Fig. 2. In the model of diesel engine, the 
turbo- charger creates a substantial coupling between 
the engine exhaust and the engine intake sides. A 
typical diesel engine model is non-linear, which is 
usually built based on equations of first physical 
principles (mass and energy) and on empirically 
adjusted equations. The nonlinear equations for each 
Diesel engine block are introduced in many studies. 
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FIGURE 1 BASIC DIESEL ENGINE CONTROL SYSTEM STRUCTURE 



<?„(Kg/h): air mass flow through the compressor. 

It (Kg/h): exhaust mass flow trough the turbine. 

qb (Kg/h): airmass flow trough the engine. 

if (Kg/h): fuel mass flow through the engine. 

1c (Kg/h): exhaust gas mass flow. 

k y (K): air temperature after the compresor. 

k 2 (K): air temperature in the intake manifold. 

£3 (K): exhaust gas temperature in front of the turbine. 

k c (K): air temperature after the compressor. 

k e (K): Exhaust gases temperature. 

Po(bar): pressure in the intake manifold. 

P\ (bar): pressure before the compressor. 

P2(bar): pressure after the turbine. 

ft (bar): pressure in the exhaust manifold. 

^(rpm.): turbocharger rotational speed. 

® e (rpm.): engine rotational speed. 

7;.(N.m): torque produced by the turbocharger. 

T e (N.m): torque produced by the engine. 

T, (N.m): resistence torque. 

"v S t (V): command signal to the VGT valve. 

\ : ratio of air and fuel flow into the cylinder. 
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FIGURE 2 THE CAUSE AND EFFECT DIAGRAM OF A TYPICAL DIESEL ENGINE 
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Wiener Model 

As shown in Fig.3, a Wiener model consists of a 
linear block followed by a static nonlinearity block, 
and it can be described as: 

A{q- l )x(k) = q- d B{q- X )u{k) (1) 
y(k) = f(x(k)) + e(k) 

Where A(q~ [ ) = l + a i q~ [ +--- + a n q~ n , 

B(q- l ) = b 0+ b iq - l +- + b m q- m . 

Where u(k) and x(k) are the input and intermediate 
variables, respectively, a, , i = 1, . . . , n and bj , j = 
1, . . . , m are the parameters of linear dynamic block. 
q~ l is the unit delay operator, and n , m are the 
orders of the linear dynamics and generally n<m . 
y(k) is the output, and /(•) is the nonlinear block of 
the Wiener model. 
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FIGURE 3 WIENER NONLINEAR DYNAMIC MODEL 

Using a polynomial function of degree p , y(k) can be 
approximately expressed as 



y(k) = f(x(k)) = Y j c,x l (k) (2) 

i=i 

Where c, , / = 1 , . . . , p is the parameter of nonlinear 
static block. I in x l (k)is the power of x(k), and p 
is the degree of the polynomial function. 

Wiener neural network for identification 

Wiener model is identified using Wiener neural 
network. Wiener neural network is designed as shown 
in Fig. 4. In the input layer of the neural network, using 
the delay operator q~ l , input signal into the system is 
expanded with time delay series nodes u(t-d) , 
u(t-d-Y),. . . , and u(t-d-m), and combined with 
feedback from the hidden layer of x(k) , the linear 
dynamic block in the Wiener model can be expressed. 
In the output layer, the Wiener model of nonlinear 
static gain block is expressed by a polynomial function 
of the hidden layer signal x(t ) . 

y(t) = fCx{t)) = Y J C l x'{k) (3) 

1=1 

Where c t , I = 1 , . . . , p is the weights, which 
associated with the parameters of nonlinear static block 
in Eq. 2. 




The hidden layer output x(f ) of the WNN can be 
expressed as 

n m 

x(t) = -^a^ik - i) + ^i j u(k -d-i) (4) 

i=l 7=1 

Where a,- , i = 1, . . . , n and bj, j = 1, . . . , m are 
the weights, which associated with the parameters of 



linear dynamic block in Eq. 1. d Is time delay of 
system? 

The task of Wiener neural network is to identify the 
parameters of a, , i = 1, . . . , n , bj, j = 1, . . . , m and 

c t ,1 = 1, . . . , p according to control input and output 
datasets. The negative gradient training algorithm can 
be used to update the weights of Wiener neural 
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network so as to obtain the identified Wiener model of 
diesel engine. The target evaluation function of the 
identifier is given as 

. n . N 

J{9j) = ™H (y(k) - m)2 = ^ Z S(fc)2 (5) 



Where y(k) and y(k) are the diesel engine output 
speed and the Wiener neural network output. N is 
the observation window length. Parameter 9 
represents the adjustable weights of a i , bj and c l . 



The gradient of the target evaluation function is 



dJ_ 
30 



2N 



dm 

80 



(6) 



According to Eqs.3, 4 and5, the partial derivatives of 
the Wiener neural network output y(k) to weights c, 

and intermediate variable x(k) can be calculated as 



dy(Jc) 
dc t 

&m_ 

dx(k) 



x' (k), I =1,2 p (7) 



= £ic,x M (Jt),/=W P (8) 



Since x(k) is also a function of a i and bj , the 
partial derivatives of intermediate variable x(k) to 
the weights a i and bj can be calculated as 



dx(k) 

dm 

dh 



-x(k — i)-^dj 

n 

u{k — d — i) — ^ a j 



dx(k - j) 
da i 

dx(k - j) 
db. 



(9) 



(10) 



From Eqs. (9) and (10), the following partial 
derivatives can be calculated as 



dy(k) = dy(k) Sxjk) ( - =12 ... n 
da i dx(k) da i 



dm _ 8m sm 

db t ~ dx{k) db. 



,i=l,2,..., m 



(11) 
(12) 



Iteration procedure of the weights can be expressed as 
is given as 

0(k + 1) = 0(k) + A0(k) = 0{k) -rj — = 6(k) + rje{k) (13) 

d3 d0 

Where ij is the training rate. 

At the beginning of each round of training, the initial 
value of the weights can be set random in range [-1, 1]. 
Iterative process can be stop, if the evaluation function 



J(8,t) is less than a set value. Eventually, the network 
weight 9 is the identification of Wiener model 
parameters. 

Bp Neural Network 

BP neural network is multi-layer perception neural 
network, which consist of three components, namely 
input layer, one or more of the hidden layer and 
output layer. Each layer has a different number of 
neurons, and the neurons of two neighboring layers 
are connected with each other. Fig.5 shows single 
hidden layer perception neural networks. 




FIGURE 5 STRUCTURE OF BP NETWORK 

In input layer, the inputs of the neurons are net) = x i , 
and outputs are O] . 

0}=x^ i=l,2,...,M ( 14 ) 

Where M is the number of the neurons in input 
layer. 

In hidden layer, the inputs of the neurons are 

M 



net 2 j (k) = ^ wfjOj , and outputs are Oj 



0)=f(net 2 m), j =\,1,...,J 



(15) 



Where wjj is the connection weight, and J is the 
number of the neurons in hidden layer. 

In Output layer, the inputs of the neurons are 



netf(k) = WpOj , and outputs are (9, 



,3 



7=1 



Of =g(netf(k)) j = 1,2,..., N 



(16) 



Where w| is the connection weight, and N is the 
number of the neurons in output layer. 

The active functions f(x) and g(x) are selected 
according to different situations. In general, the 
activation functions are both Sigmoid-function. 



fix)-- 



(17) 
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8(x)=- 



e +e 



(18) 



In Eq.20, 



PID Controller Based on Improved BP Neural 
Network and Wiener Neural Network 

As shown in Fig.6, an adaptive controller is proposed 
to control speed of diesel engine, which includes a 
Wiener model identifier and a PB-PID controller. The 
PID controller based on BP neural network (PB-PID) 
consists of a PID controller and a BP neural network. 
The PID parameters are the three outputs of BP neural 
network. According to the running status of control 
object, PID parameters can be adjusted by neural 
network. 

PID control is often used to control diesel engine speed. 
It is a negative feedback control, and digital 
incremental form of traditional PID controller can be 
described as 

Au{k) = K p (e(k) - e(k - 1) + K ,e{k) + K D (e{k) - 2e(k - 1) + e{k - 2)) 
Where K p is the proportional coefficient. Kj is 
the integration time ,and K D is the differential 
time. 
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FIGURE 6 PID CONTROLLER BASED ON BP NN AND WIENER 

NN 

Considering the relative speed error at different 
working speeds, relative error is used in target 
evaluation function of the BP neural network. The 
target evaluation function of BP-PID is defined as: 



E(k): 



(r m (k)-you,(k)r 
2- !*,(*) 



(19) 



The purpose of learning algorithm of BP network is to 
minimize E(k) . The gradient descent algorithm is used 
to modify weights of BP network. The connection 
weight of output layer is adjusted as follow. 

8E(k) 



Aw? ] {k) = -r 1 



+ 8Aw?Ak-\) (20) 
dwl(k) 



Where rj # is learning rate, and # d < 
coefficient. 



is inertia 



8E(k) _ dE{k) dy out {k) 8u(k) 80? (k) dnetf(k) 



dw?j(k) dy out (k) du(k) 80? (k) dnetf(k) dwfj(k) 

= _£W.^). W ■ 3 2 (21) 
r m du(k) 80? (k) ' 

If O? (k) = K p , 0\ (k) = K, , Ol (k) = K D , 
8u(k) 



80^(k) 
8u(k) 



-e{k)-e{k-\) 



8u(k) 
dOjik) 



= e(k) 
80 3 2 {k) 

-e(k)-2e(k-\) + e(k-2) 



Assuming, 



S?=e(k). 8y °»< (k) -^ 
8u(k) 80? (k) 

Where g Q = g(x)(l- g(x)) . 

Eq. 20 can be expressed as follow. 



Aw?. (k) = ? 7 ^-Oj(k) + 8 Awl {k ~ J ) 



■g\net?(k)) (22) 



For the same method, 



Aw 2 H (k) = ri-!-net) (k) + 8Aw 2 H (k - 1) 



(23) 



(24) 



j wh ere S? = ^ w| (k) ■ f ' (net] (k)) , and /'(•) = (1 - / 2 (*)) / 2 . 



8y out (k)l 8u(k) denotes the sensitivity of the diesel 
engine output with respect to its input. Because 
8y out {k)l 8u(k) is unknown, it replaced by 
sgn [8y(k) 1 8u(k)] usually. The sensitivity cannot be 
calculated directly from the output of the nonlinear 
system since the precise mathematical model is usually 
unknown in a diesel engine control system. In general, 
since the Wiener neural network was trained, we can 
use Wiener model to displace the diesel engine control 
model. Namely, y out (k) = y out (k) . Therefore, the 
system sensitivity can be calculated as 

dy u,( k ) g dy „t( fc ) = dy out (k) 8x(k) 
8u(k) 8u(k) 8x(k) 8u(k) 

Experiment of the Adaptive Controller 

In order to validate the effectiveness of the adaptive 
controller, an experiment of control speed of a diesel 
engine was done. The experiment is based on the 
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dynamic simulation model of a D6114 diesel engine 
produced by Shanghai, which was modeled with 
exhaust-filled method. The model is calibrated with 
actual steady state datum. In the adaptive controller, 
Wiener neural network to identify the diesel engine 
model. The initial parameters in Wiener neural 
network are set to be y(k) = 0, x(k)= 0, dy(k)/dc,=0, 
dx(k)/dai=0, and dx(k)/dbt=0 when k < 0. The initial 
values of the weights a ir b jr andc, and are all selected 

within [-1, 1] randomly. The training rate is selected as 
77=0.5. Choosing n =3, m=2, p=3, J < 0.01, and the 
maximum number of iterations be 200. 



In Fig.7, the closed-loop responses are shown for the 
proposed adaptive BP-PID and PID controllers. The 
diesel engine rotational speed command varied 
alternately from 1100 (RPM) to 1500 (RPM). The figures 
illustrate that the diesel engine controlled by the 
BP-PID controller has a better tracking performance 
than using the PID controller. Controlled by the 
BP-PID controller, the speed fluctuations rate is less 
than 0.53%, and stability time is not more than 4 
seconds; Controlled by the PID controller, the speed 
fluctuations rate is about 0.8%, and stability time is 
more than 5 seconds. 
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FIGURE 7 RESPONSE CURVE OF THE PB-PID AND PID CONTROL SYSTEMS 



Conclusion 

In this study, an improved BP neural network is used 
to tuning PID parameters of a diesel engine speed 
controller. Considering the speed error at different 
working speeds, relative error is used in target 
evaluation function of BP neural network. Moreover, a 
Wiener neural network is used to identify the model 
of the diesel engine in order to improve the rate of 
convergence of BP neural network. The results of the 
simulation experiments show that the PID controller 
based on the improved BP neural network has faster 
speed, stronger adaptability. 
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